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We introduce density dependence of the cell size in cellular-automaton models for traffic flow, 
which allows a more precise correspondence between real- world phenomena and what observed 
in simulation. Also, we give an explicit calibration of the particle density particularly for the 
asymmetric simple exclusion process with some update rules. We thus find that the present 
method is valid in that it reproduces a realistic flow-density diagram. 
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exact solution 



In the last several decades, simulation research has 
firmly established its place in physics. Especially, cellular 
automaton (CA) modeling plays a central role in hydro- 
dynamics, statistical mechanics and general phenomcno- 
logical studies including complex systems. 1-3 ) CA mod- 
els are, naturally, suitable for simulation on a computer, 
but there are also actual discrete systems like CA, e.g., 
molecular motors on a microtubule. 4 ' In such systems, 
particles have a finite volume, and in principle, both col- 
lision and overtaking are forbidden. As a result, we often 
observe some particles stemming the flow by themselves. 
The above situation can be precisely mimicked with the 
so-called exclusion rule in CA models. 

Introduction of CA models produced a breakthrough 
in research on traffic flow, 5-12 ) and now CA models 
are widely applied in real- world technologies, such as in 
real-time traffic forecast. One of the most basic models 
for traffic flow is the asymmetric simple exclusion pro- 
cess 13 ' 14 ) (ASEP), a CA model incorporating random- 
ness into an elementary CA rule 184; 1 ) particles in the 
ASEP hop to the next site with a given probability if it 
is not occupied. We should remark here that there are 
some possible choices to update cells, 15 ' 16 ) e.g., random 
sequential updating and shuffled dynamics, 17 ) whereas 
all sites are synchronously updated in a CA. 1 ) However, 
we do not adhere to only the parallel update case in the 
following discussion. 

The validation of traffic-flow models is carried out 
by comparing the simulation result to observational 
data. 7 ' 9 ' 10,18-21 ) Normal traffic data contains the den- 
sity, velocity and flow, and one usually plots the velocity 
against the flow, the velocity against the density and 
the flow against the density. In particular, the signifi- 
cant properties of a traffic flow are summed up in the 
flow-density diagram. 22 ) Some simple models such as the 
ASEP are known to be exactly solvable, 15 ' 16 ' 23 ' 24 ) which 
means that one can obtain an exact flow-density dia- 
gram. Note that the flow-density diagram for the ASEP 
has the so-called hole-particle symmetry: i.e., inverting 
all the cells in concert, one sees that the diagram is sym- 
metric in the density. In what follows, we discuss a cal- 
ibration of CA models for traffic flow focusing on these 
exactly solvable ones. 
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PCA = -r (0 < Pca < 1) 



(2) 



In this work, we simply consider v vehicles of the same 
size c on a circuit of length I. The real- world density of 
vehicles, which is conserved in simulations, is accordingly 
defined as 

vc . 

Prw = — (0 < Prw < 1)- (1) 

Meanwhile, in the corresponding CA model the particle 
density is defined as 

N 
T 

where N is the number of occupied cells, and L is the 
total number of cells. Let d be the cell size. Since v = N 
and I = Ld, we then have 

Pca _ d 
Prw c 

In this letter, we shed light on the following condition 
with respect to vehicle size c and cell size d: 

c<d<2c. (4) 

This is justified because if d < c then vehicles could not 
be contained in each single cell, or if 2c < d then one 
should reduce the cell size d to half. The CA models for 
traffic flow admit of optimizing the cell size according 
to the particle density. We therefore assume that the cell 
size d is determined depending on the density p RW as d = 
c/(pkw)j where f(p) is a function such that 1 < /(p) < 2 
due to eq. (4). Accordingly, eq. (3) yields 



(3) 



Pca — Prw/(Prw)- 



(5) 



Note that since < p CA < 1 by definition, we have to 
refine the condition on /(p) to 



1 < /(p) < min{2, l/p}. 



(6) 



One may use a larger cell size while the density is low, 
but a smaller one while it is high. Thus, wc define /(p) 
as a non-increasing function in p. For example, a simple 
choice, f(p)= 2 — p, leads to 



Pca — Prw (2 — Prw)- 



(7) 



As mentioned above, the ASEP is exactly solvable. The 
flow Qrandom is described as a function of the density p CA : 



i(Pca) = PPca(1 - Pca), 



(8) 
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Fig. 1. The ASEP with random sequential updating and hop 
probability p = 1: (a)The flow-density diagram and (b)The 
velocity-density diagram; (solid line) calibrated and (dashed 
line)non-calibrated. 

where p is the hop probability. Note that eq. (8) is for 
the random sequential updating. Then, using eq. (7), the 
flow-density diagram, eq. (8), is calibrated as 

Qrandom(PRw) = PPrw{2 — /?Rw)(l — Prw) 2 - (9) 

(Note that we often denote different functions by the 
same symbol for the sake of simplicity.) In Fig. 1(a), we 
show the graphs of eqs. (8) and (9), finding that the 
calibration (7) helps to reproduce a more realistic flow- 
density diagram; in particular, the maximum-flow den- 
sity moves from 1/2 to 1 - l/y/2 (~ 0.29). 

Also, we consider the average velocity of vehicles as 
a function of the density. Following the assumption, 22 ) 
Q = pv, with respect to flow Q, density p, and velocity 
v, we have 

Wrandom(PCA) =P(l - PCa) , ^ 
^random (pKXv) =J*(2 — Prw)(1 — Pkw) • 

Figure 1(b) shows the graphs of eq. (10). We see density 
dependence of the mean hop probability in the ASEP. 
Note that after the calibration the average velocity is no 
longer equivalent to the mean hop probability and can 
take values larger than 1. 

We add another example: the ASEP with the parallel 
updating which includes the rule- 184 CA as a special 
(p = 1). The flow-density diagram for the parallel update 
case, which is described as 23 ) 

Qparallel(PcA) = 5 t 1 _ V 1 ~ 4 PPca(1 ~ Pca)] (H) 
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Fig. 2. The ASEP with the parallel updating and hop probability 
p = 0.9, and those of the rule-184 CA (p = 1): (a)The flow- 
density diagram and (b)The velocity-density diagram; (black 
solid line)ASEP calibrated, (gray solid line)Rule-184 CA cali- 
brated, (black dashed line)ASEP non-calibrated, (gray dashed 
line)Rule-184 CA non-calibrated. 

is calibrated, using eq. (7), to be 
„ , ^_ 1 - V 1 ~ 4 PPrw(2 - p RW )(l - Prw) 2 

V parallel (,/fRW J — ~ ' 

(12) 

In Fig. 2(a), we show the graph of eq. (12). The 
maximum-flow density calibrated is equal to 1 — l/v2, 
and is independent of p. Moreover, it is identical to that 
for the random sequential update case. In the special case 
of the rule-184 CA, we have 

Qparallel(/3CA) = min{/5 C A, 1 ~ PC\} , 

(13) 

Qparallcl(/3Rw) =mm{/9 HW (2 - jO R w)) (1 - Pltw) 2 } ■ 

In this case (p = 1 and the parallel updating), we see a 
first-order phase transition at the maximum-flow density, 
which distinguishes the model from the other ones. The 
phase transition implies that a traffic jam never occurs 
until the density exceeds the maximum-flow density. (By 
contrast, we find a traffic jam occurring at rather lower 
densities in the other cases.) 

As well as in the previous case, the average velocity is 
obtained from the flow via v = Q/p. In particular, from 
eq. (13), those for the rule-184 CA are given by 

«parallcl(/0 C A) = min{ 1 , -1 + 1//? 0A }, 

(14) 

Uparallel(PRw) =mm{2 - p RW , (1 - /) RW ) 2 /ftw}- 

Figure 2(b) shows the graphs of eq. (14), which presents 
how the proposed calibration works. Once the cell size is 
fixed, dynamics of the real-world vehicles is transformed 



J. Phys. Soc. Jpn. 



Letter 



Author Name 3 



into a coarse-grained motion in the CA model: Particles 
therein hop at most one site, with a constant probability 
and the exclusive interaction. The calibration then helps 
one restore the original dynamics. Compare Figs. 1 and 2 
with the simulated/observed diagrams given in refs. 7, 8, 
21, 23, and 24 to see the effect of the present calibration. 

In this letter, we proposed a calibration of the particle 
density in CA models for traffic flow. One of our motiva- 
tions in the present work is to reduce the cost of large- 
scale traffic simulations such as city traffic, using some 
simple model which reproduces a realistic traffic flow. 
Our central idea is that the cell size, which is usually set 
to the vehicle size, may also depend on other parameters, 
e.g., the number of vehicles conserved in each simulation. 
It is not difficult to apply this idea for general CA mod- 
els such as two-dimensional CA models for pedestrian 
flow and computation techniques (e.g., the Monte-Carlo 
method). We need however to modify the models in or- 
der to reproduce detailed structures of traffic flow, such 
as a metastable branch appearing in the observed flow- 
density diagram. Nevertheless, as far as traffic flow is re- 
garded as a fluid, i.e., on a macroscopic level, these simple 
CA models share a lot with a real-world traffic flow. We 
conclude, more precisely, that with a calibration of the 
particle density, the exclusion rule in CA models shows, 
beyond the hole-particle symmetry, an asymmetric flow- 
density relation peculiar and fundamental to traffic flow. 
Finally, it is fair to comment that the calibration method 
proposed works reliably only for homogeneous systems 
such as a periodic lattice. However, we believe that we 
modify the present method so that one can apply it to 
the ASEP with an open boundary condition. It will be 
reported in our subsequent publications. 
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